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singularities in i* 1 (G)-gravity and J~(R, GQ-gravity by taking into account higher-order curvature 
corrections. The example of non-singular realistic modified Gauss-Bonnet gravity is presented. 
It turns out that adding such non-singular modified gravity to singular Dark Energy makes the 
combined theory to be non-singular one as well. 
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I. INTRODUCTION 



Recent observations have implied that the current expansion of the universe is acceler- 
ating |2|. There exist two broad categories to explain this phenomena UQ. One is the 
introduction of "dark energy" in the framework of general relativity. The other is the investi- 
gation of a modified gravitational theory, e.g., f{R) -gravity, in which the action is described 
by the Ricci scalar R plus an arbitrary function f(R) of R (for reviews, see 6l-ll0|). 

It is known that accelerating Friedmann-Robertson- Walker (FRW) universe is described 
by cosmological constant /quintessence/phantom Dark Energy. In principle, Dark Energy 
(DE) could be described by scalar field theories, fluid, modified gravity, etc. It is quite 
well-known that any of such DE models may be represented as the effective fluid with 
corresponding characteristics. At the late-time accelerating stage of the FRW universe, if 
the ratio of the effective pressure to the effective energy density of the universe, i.e., the 
effective equation of state (EoS) w c g = p e s/p e s, is larger than —1, it is the quintessence 



15] (non-phantom) phase. On the other hand, if w e g is less than —1, it is the phantom 



phase 



16] while effective cosmological constant appears as DE when w e ff 



Note 



that (non-transient) phantom phase evolution usually ends up in Type I (Big Rip) future 
singularity. It is remarkable that many of the effective quintessence/phantom DEs may 
bring the future universe evolution to finite-time singularity. The classification of such 
finite-time future singularities has been made in Ref. [17]. Some of these four types future 
singularities are softer than other, for instance, not all characteristic quantities (scale factor, 
effective pressure and energy-density) diverge in rip time. There is not any qualitative 
difference between convenient DEs and modified gravity in this respect. For instance, the 
convenient parameter-dependent DE models may show all four possible types of finite-time 



future singularity as demonstrated in Refs. 



in Refs. 120 
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19] . On the same time it was demonstrated 



2l] that /(-R)-gravity DE model may also bring the universe evolution to all four 



possible future sing ularities (for first example of Big Rip (Type I) singularity in modified 



gravity, see 



22 



24j). Furthermore, it is interesting that f(R) modified gravity may also 
provide the universal scenario to cure the finite-time future singularity by adding, say, R 2 - 
term 20|, |2l|, |23|, |25| or non-singular viable /(i?) -gravity 26] (for related discussion of Type 
II future singularity in particular /(i?)-gravity and its curing by i? 2 -term, see Refs. 0-0 )• 



It is clear that singular dark energy may lead to various instabilities in the current universe 
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cosmology, including black holes and stellar astrophysics. In this respect it is very important 
to list the singular dark energy models as well as try to indicate the physical consequences 
of possible future singularity. Moreover, it is desirable to construct the universal scenario to 
cure such singularities. The interesting class of modified gravity models which may easily 
produce the late-time acceleration epoch is string-inspired modified Gauss-Bonnet gravity, 



so-called F(G)-gravity 33M36], where F(G) is an arbitrary function of the Gauss-Bonnet 



invariant G = R 2 - AR^RT + R^R^ (R^ and R^ are the Ricci tensor and the 
Riemann tensor, respectively). It is known that such class of models may also lead to 



finite-time future singularity 



2i 



2l| . we explore the F(R, G)- 



In the present paper, as a generalized investigation of Ref. 
gravity models with realizing the finite-time future singularities by using the reconstruction 

37| . where F(R, G) is an arbitrary function of R and G. 



method of modified gravity 
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The F(R, G)-gravity is a gravitational theory in a more general class of modified gravity 
and includes /(i?)-gravity and F(G)-gravity. As an explicit example of F(R, G)-gravity, we 
also investigate F(G)-gravity and reconstruct the F(G)-gravity models in which finite-time 
future singularities could appear in great detail. It is shown that all four types of finite- 
time future singularity may occur in such modified gravity. In addition, we examine the 
possibility of the finite-time future singularities in F(G)-gravity and F(R, G)-gravity being 
cured under higher-order curvature corrections. The explicitly non-singular modified Gauss- 
Bonnet models is proposed and it is shown that the finite-time future singularities may be 
easily protected combining a singular theory with the non-singular one. This suggests the 
universal scenario to cure the finite-time future singularity in the same line as it was proposed 
in Ref. [26 1. 



The paper is organized as follows. In Sec. II, we briefly review the model of F(R, G)- 
gravity and write down the gravitational field equations. In addition, we classify the four 
types of the finite-time future singularities. In Sec. Ill, as a first step, we investigate F(G)- 
gravity and reconstruct the F(G)-gravity models where finite-time future singularities may 
occur. We also examine the finite-time future singularities in realistic models of F(G)- 
gravity. Next, in Sec. IV we study the general F(R, G)-gravity models where the finite-time 
future singularities occur. Moreover, we explore the finite-time future singularities in a 
realistic model of F(R, G)-gravity. In Sec. V, we discuss a possible way to resolve the 
finite-time future singularities in F(G)-gravity and F(R, G)-gravity by taking into account 
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higher-order curvature corrections. The non-singular theories are proposed. It is shown 
that the addition of such non-singular effective dark energy to the singular one may cure 
the singularity of the combined theory. Hence, modified Gauss-Bonnet gravity may appear 
as the effective universal regulator of finite-time future singularity not only for singular 
alternative gravity but also for convenient singular DE. Finally, conclusions are given in Sec. 
VI. The finite-time future singularities in a simple model of /(i?)-gravity are also examined 
in Appendix A. Furthermore, a further argument on the asymptotic behavior of singular 
models is presented in Appendix B. 

We use units of k-B — c — H — 1 and denote the gravitational constant 8ttGn by k 2 = 
8n / Mpi 2 with the Planck mass of Mp\ = G^ 1 ^ 2 = 1.2 x 10 19 GeV. A note on notation is that 
throughout the present paper, a, 7, z, n, m, 5 and ( are constants unless we mention some 
conditions in regard to these expressions. 



II. T(R, GO-GRAVITY 

In this section, we briefly review J-"(-R, G)-gravity and derive the gravitational field equa- 
tions. Moreover, we classify the finite-time future singularities into four types following 
ref. |l7(|. 



A. The Model 



The action of J-"(-R, G)-gravity is given by 

~F(R, G) 



S 



2k 2 



matter 



(2.1) 



where g is the determinant of the metric tensor g pv and £ ma tter is the matter Lagrangian. 
From the action in Eq. (12. ip . the gravitational field equation is derived as 

- \Rg»)\ = ^XT tter) + (J 7 - FrR) + V^V^ - g^UT' R 
+ (—2RR pu + AR pp R v p — 2R p paT R upaT + Ag af 'g^ c 7 R^avpRpa)^ 
+2 (V M V„J^) R - 2 9lxu (UF'c) R + 4 (DT' G ) - 4 (V P V^) R v p 
-4 (y p V v F' G ) V + A 9flP (V p V a F' G ) BT - 4 (V p V a F G ) g ap g^R, a ^ , (2.2) 
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where we have used the following expressions: 

_ dF{R, G) _ dJ r (R, G) 

>i? " dR ' dG • {2 - 6) 

Here, V M is the covariant derivative operator associated with g^, □ = g^V^Vu is the 

covariant d'Alembertian for a scalar field, and j^™ atter ) — fc^g (p ; p ; p ; p) i s the contribution 

to the energy-momentum tensor from all ordinary matters with p and p being the energy 

density and pressure of all ordinary matters, respectively. 

The most general flat FRW space-time is described by the metric 

ds 2 = -N 2 (t)dt 2 + a 2 (t)dx 2 , (2.4) 

where a(t) is the scale factor and N(t) is an arbitrary function of t. In what follows, we take 
N(t) = 1. 

In the FRW background, from (p, v) = (0, 0) and the trace part of (/i, v) = = 
1, • • • ,3) components in Eq. (12. 2p . we obtain the gravitational field equations: 



PcS = ^H 2 , PcS = -^(2H + 3H 2 ) , (2.5) 



where p eS and p e s are the effective energy density and pressure of the universe, respectively, 
and these are defined as 



Peff = -^r 

•'R 



PeS = JF< i P+ 7^2 

R 



{ P + 2^ ~' F) ~ QHPr + GT ' G ~ 2AH " P <] } ' (Z6) 

{ P + 2^ [~ ^' rR 4HPr + 2Pr ~ GTg + 16H (H + # 2 ) P g 



2-cl 

G 



(2.7) 



Here, H = a(t)/a(t) is the Hubble parameter and the dot denotes the time derivative of d/ dt. 
For general relativity with T{R, G) = R, p c g = p and p e g = p and therefore Eqs. ( 12. 6ft and 
( 12. 7p are the FRW equations. Consequently, Eqs. (12. 6p and (12. 7p imply that the contribution 
of modified gravity can formally be included in the effective energy density and pressure of 
the universe. 



B. Four types of the finite-time future singularities 



We consider the case in which the Hubble parameter is expressed as 



where h, to and Hq are positive constants, (3 is a constant, and t < to- We can see that if 
/3 > 0, H becomes singular in the limit t — > to. Hence, to is the time when a singularity 
appears. On the other hand, if /3 < 0, even for non-integer values of /3 some derivative of 
H and therefore the curvature becomes singular |21j. We assume /3 ^ because /3 = 
corresponds to de Sitter space, which has no singularity. 

The finite-time future singularities can be classified in the following way Jl7| : 

• Type I (Big Rip): for t — > t , a(t) — > oo, p e g — > oo and \p c s\ — > oo. The case in which 
Peg and p e ff are finite at to is also included. It corresponds to = 1 and (3 > 1. 

• Type II (sudden 38]): for t — > t , a(t) — > Oq, p e fr - > Po and |p eff | — oo. It corresponds 
to -1 < < 1. 

• Type III: for £ — > t , a(t) — > cto, p c ff —> oo and \p c s\ — )■ oo. It corresponds to < /3 < 1. 

• Type IV: for t — > to, a(t) — > ao, p c s 0, \p c s\ — > and higher derivatives of H diverge. 
The case in which p and/or p tend to finite values is also included. It corresponds to 
f3 < —1 but (3 is not any integer number. 

Here, a (^ 0) and po are constants. We note that in the present paper, we call singularities 
for f3 — 1 and those for (3 > 1 as the "Big Rip" singularities and the "Type I" singularities, 
respectively. 



III. F(G)-GRAVITY 



33 



In this section, as an explicit example of J-~(R } G)-gravity, we first study F(G)-gravity 
36]. We reconstruct the -F(G)-gravity models where finite-time future singularities may 
occur. In addition, we explore the finite-time future singularities in realistic models of 
F(G)-gravity. 



A. The Model 



The action of F(G)-gravity is given by 33| 

1 



S = Id Xyf^g 



2k 



_ (R + F(G)) + C 

matter 



(3.1) 
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which corresponds to the action in Eq. (12. ip with J-~(R, G) = R + F(G). 

In the FRW background in Eq. (12.41) with N(t) = 1, it follows from the action in Eq. (13. ip 
that the equations of motion (EOM) for F(G)-gravity are given by 35| 



2AH 3 F'(G) + 6H 2 + F(G) - GF\G) = 2k 2 p , (3.2) 
8H 2 F\G) + 1QHF'(G) (h + H 2 ) + (aH + 6H 2 ) + F(G) - GF'(G) = -2n 2 p , (3.3) 



where the prime denotes differentiation with respect to G. Moreover, we have 

R = 6 (2H 2 + H 
G = 2AH 2 (H 2 + H 



(3.4) 
(3.5) 



In this case, p e s and p e g in the FRW equations (12.51) take the form 
1 



PcS 



2k 2 
+P- 



—F(G) + 2AH 2 (H 2 + H) F'(G)-2A 2 H 4 ( 2H 2 + H H + 4H 2 H) F"(G) 

(3.6) 

QH 3 + 8HHH + 2AH 2 H 2 + QH 3 H 

(3.7) 



p cS = ^{ F ( G ) - 24H \ H +H)F'(G) + (24)8tf 
+8H 4 H + H 2 H 



F"{G) + 8{2A) 2 H i 2H 2 + HH + AH H F"'(G) \ + p 



where we have used Eq. (13. 51) . 

We assume that the matter has a constant equation of state (EoS) parameter w = p/p. 
By combining the two equations in Eq. (12.51) . we obtain 



G(H,H. 



2H + 3(1 + w)H 2 



where 



g(H,H...) = p cH - wp, 



cff 



(3.8) 



(3.9) 



When a cosmology is given by H = H(t), the right-hand side of Eq. (13. 8p is described by 
a function of t. If the function Q(H,H...) in Eq. (13. 9p . which is the combination of H, H, 
H and the higher derivatives of H, reproduce the above function of t, this cosmology could 
be realized. Hence, the function Q(H,H...) can be used to judge whether the particular 
cosmology could be realized or not [2l|. The form of Q(H, H...) is determined by the gravi- 
tational theory which one considers. In the case of F(G)-gravity, by substituting Eqs. f ]3.6f) 
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and (E2D into Eq. flU), we find 



G(H, H...) = ^2 1(1 + w)F(G) - 24(1 + w)# 2 (# 2 + #) F'{G) + 8(24)# 2 

r2 rr2 



QH 3 



+8HHH + 6(4 + w)H 2 H 2 + 3(2 + w).FP# + 4(2 + 3w)H A H + H H 



+8(24) 2 # 4 ( 2if 2 + HE + 4H' 2 H ) F"'(G) } . 



2 u 



F"(G) 
(3.10) 



B. Finite-time future singularities in F(G)-gravity 



We investigate the F(G)-gravity models in which the finite-time future singularities could 
occur, when the form of H is taken as Eq. (12.81). T o find such F(G)-gravity models, we use 
the reconstruction method of modified gravity (21I . By using proper functions P(t) and 
Q(t) of a scalar field t which we identify with the cosmic time, the action in Eq. (13.11) can 
be rewritten to 



S 



d A Xy/—g 



^(R + P(t)G + Q(t)) + £ mat tcr 



The variation with respect to t yields 



dP(t\ G dm = 



(3.11) 



(3.12) 



dt dt 

from which we can find t = t(G). By substituting t = t(G) into Eq. (|3.11j) . we find the 
action in terms of F(G) 

F(G) = P{t)G + Q(t) . (3.13) 

We describe the scale factor as 



a(t) = aexp (g(t)) 



(3.14) 



where a is a constant and g(t) is a proper function. By using Eqs. (13. 2J) , Eq. (13.31) . (I3.14p . 
the matter conservation law p + 3H(p + p) = and then neglecting the contribution from 
matter, we get the differential equation 



d 



dP(t)\ 



dP(t) 
dt 



+ g(t) = o. 



By using the first EOM for F(G)-gravity in Eq. (|3.2|) . Q(t) is given by 

Q(t) = -24«7 3 (t)^ - 6g\t) . 



(3.15) 



(3.16) 
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1. Big Rip singularity 



First, we examine the Big Rip singularity. If J3 = 1 in Eq. (12.81) with Hq = 0, H and G 
are given by 

h 



H 

G 



(to - 1) 
24h 3 



(to - ty 

The most general solution of Eq. (13. 15|) is given by 



1 + h). 



P(t) 



1 , , (to~tf~ h 

(2t - t)t + — '- + c 2 



Ah(h-iy~~" * 3-h 

where c\ and c 2 are constants. From Eq. ( I3.16p . we get 



(3.17) 
(3.18) 

(3.19) 



Q(t) 



6h 2 



24h 3 



(to-t) 
2h(h-l) 



- Ci(to - *) 



2-h 



(t - t)> (t - ty 

Furthermore, from Eq. (13. 12[) we obtain 

~ 24(h 3 + h 4 ^ 1/4 
G 



(3.20) 



(3.21) 



which is consistent with Eq. (13.181) . By solving Eq. (I3.13p . we find the most general form of 
F(G) which realizes the Big Rip singularity 



J6h 3 (l + h) r- „H+X 

F(G) = v , lA , N vG + c\G 4 + c 2 G. 



(3.22) 



This is an exact solution of Eq. ( 13. 8 j) in the case of Eq. ( I3.17p . In general, if for large values 
of G, -F(G) ~ aG 1 / 2 , where a(^ 0) is a constant, the Big Rip singularity could appear for 
any value of h ^ 1. In the case of h = 1, the solution of G(H, H...) is zero for F(G) = aG 1 ^ 2 . 
Note that aG^ 1+h ^ 4 is an invariant with respect to the Big Rip solution. 
In the case of h — 1, it is possible to find another exact solution for P(t) 



P(t) = a(t -t)nn[ 1 (t -t) z ] 



(3.23) 



where 7(> 0) is a positive constant and q and z are constants. The equation ( 13. 15ft is 
satisfied for the case of Eq. (I3.17P if g = 3 — h = 2 (and therefore h — 1) and za = —1/4. 
From Eq. (I3.16p . we have 

12 



Q(t) 



(to-ty 



In b(t -t)) 



(3.24) 



The form of F(G) is given by 



F{G) = —VG\n(jG). (3.25) 



This is another exact solution of Eq. ( 13.81) for H = l/(to — t). In general, if for large values 
of G, F(G) ~ ay/Gln^G) with a > and 7 > 0, the Big Rip singularity could appear. 
The same result is found for F(G) ~ ay/G In^G 2 + Gq) with a > 0, 7 > and z > 0, where 
Go is a constant. 

2. Other types of singularities 

Next, we investigate the other types of singularities. If /3 7^ 1, Eq. (12. 8p with Hq = 
implies that the scale factor a(t) behaves as 

'h(t -ty~e- 



a(t) = exp ^ 
We consider the case in which H and G are given by 



(3.26) 



ff= (vV (3 - 27) 

24/? 4 

G ~<^r' (3 - 28) 



A solution of Eq. ( 13. 151) in the limit t — > to is given by 



P W - TT^I ( 3 - 29 ) 
(*o - t) 

with ^ = —2/3 and a = — l/4/i 2 . The form of F(G) is expressed as 

F(G) = -12y^. (3.30) 

Hence, if for large values of G, F(G) ~ —ay/G with a > 0, a Type I singularity could 
appear. 

When /3 < 1, the forms of H and G are given by 

24fr 3 /3 

G ~ (t - t)3/m ' (3 - 32) 
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An asymptotic solution of Eq. (I3.15P in the limit t — > to is given by 

P(t) - ^ (3.33) 
with z = — (1 + 0) and a = l/2/i(l + 0). The form of F(G) becomes 

a h 2 / i r i \ 2/3/(3/3+1) 

^(G0 = 7^TT(3/? + l) UlT^T • (3-34) 



+ 7 V 24/i3 

Hence, if for large values of G, F(G) has the form 

F(G)~a\G\\ 7 = ^-7, (3.35) 

with a > and < 7 < 1/2, we find < < 1 and a Type III singularity could emerge. 

If for G — > —00, ^(G) has the form in Eq. (13. 35ft with a > and —00 < 7 < 0, we find 
— 1/3 < < and a Type II (sudden) singularity could appear. Moreover, if for G — > 0~, 
F(G) has the form in Eq. (I3.35P with a < and 1 < 7 < 00, we obtain — 1 < < — 1/3 and 
a Type II singularity could occur. 

If for G (T, F(G) has the form in Eq. fl3~35|) with a > and 2/3 < 7 < 1, we 
obtain —00 < < — 1 and a Type IV singularity could appear. We also require that 
7 7^ 2n/(3n — 1), where n is a natural number. 

We can generate all the possible Type II singularities as shown above except in the case 
= —1/3, i.e., H = h/(to — t) 1 ^ 3 . In this case, we have the following form of G: 

G = 24h 3 + 24/i 4 (t - t) 4/3 < . (3.36) 

To find t in terms of G, we must consider the whole expression of G by taking into account 
also the low term of (t — t). We obtain 

F(G) ~ -^^G{G + Sh 3 ) 1 / 2 + 2=(G + 8hY 2 , (3.37) 

which satisfies Eq. (13. 8p in the limit t — > to- As a consequence, the specific model F(G) = 
<j\G(G + C3) 1 / 2 + cr 2 (G + C3) 1 / 2 , where a±, 02 and C3 are positive constants, can generate a 
Type II singularity. 
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C. Realistic models of F(G)-gravity 



Here, we study the realistic models of F(G)-gravity, which reproduce the current accel- 
eration, namely [36| 

= < 339 > 

F 3 (G) = a 3 G n (l + b 3 G m ), (3.40) 

where ai, 02, b%, 62, 03, &3, ?i, ^ and m are constants. In the following, we always assume 
n > 0. For the model (I3.39p . Types I, II and III singularities may be present. In fact, for 
N = 1/2, one could have Big Rip singularities, since in this case, in the limit large G, Eq. 
(I3.39P gives aG 1 ^ 2 . Thus, as discussed in Subsection III. B, one has a Big Rip singularity. 
Moreover, again with iV = 1/2, if 01/02 < 0, Eq. ( 13.391) for large value of G, leads to 
— olG 1 ! 2 with a > and thus Type I singularity could appear. If n and N are integers and 
n + N > 0, for large and negative value of G, FziG) ~ axja^G . As a result, a Type II 
singularity could appear, when — n < N < 0, N even and a\/a 2 > or N odd and ai/a 2 < 
(see Eq. (I3.35P and the related discussion). If < iV < 1/2 and 01/02 > 0, we have the 
Type III singularity (see Eq. (I3.35P ). When G — > 0~, we do not recover any example of 
singularity of the preceding subsection. 

If there exists any singularity solution, it must be consistent with Eq. f)3.8p . The behavior 
of Eq. (13.1 Op takes two asymptotic forms which depend on the parameter of as follows: 

• Case of (3 > 1: In the limit t — > to, we find 

g(H, *...)- aF(G) + j^F> ( G) + ^^^"(0) + j^HG) , 

(3.41) 

where 5 and £ are constants. To realize a singularity, from Eq. (I3.8P we must have 

• 3(1 + w)h 2 

g(H,H...)~--^-—^. (3.42) 

Hence, if for G ~ 24o 4 /(t - t) 4/s with f3 > 1, the highest term of Eq. (13^B is 
proportional to l/(to — t) 2/3 , it is possible to have a Type I singularity. This condition 
is necessary and not sufficient. Another very important condition that must be satisfied 
is the concordance of the signs in Eq. (13. 42 p . which depends on the parameters of the 
model. 
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Case of /3 < 1: In the limit t — > to we obtain 

(to -V +lF ' (G) + '(*n (/..-/) 



if...) ~ aF(G) + ^ _\ W+1 F'{G) + f+ __ 6 A50+3 F"(G) + f+ _ _is^ F '"( G ) ■ 



(3.43) 

To realize a singularity, from Eq. (13. 8p we must have 

(3 ' 44) 

Thus, if for G ~ 24h 3 f3/(t - t) 3/3+1 with < 1, the highest term of Eq. (13T431 is 
proportional to l/(to — t)^ +1 , it is possible to have a Type II, III or IV singularity. 
Also this condition is necessary and not sufficient. 

We see that the model in Eq. (13.391) with n > and N > is not able to realize a Type 
IV singularity because for /3 < — 1 the right-hand side of Eq. (I3.44p tends to zero and the 
left-hand side of Eq. (13.44)) tends to a constant (F^iG) ~ bx/b^). Nevertheless, it is possible 
to have a Type II singularity for < /3 < — 1/3. If n > and N > 0, we get 



f a (G)~£, F&G) ~ -n-^G^\ F%{G) ~ -n-^r(n - l)G n ~' 2 



'2 



F%{G) ~ - l)(n - 2)G n - 3 . (3.45) 



y 2 

It can be shown that, under the requirement n > 1 (the relation between n and /3 is 
n = 2/3/(3/3 + 1)), the asymptotic behavior of Eq. (g3SD when G ~ 24/i 3 /3/(t - t) 3/3+1 is 
proportional to l/(io — t) /3+1 and therefore it is possible to realize the Type II singularity. 

• For N —1 and n = 2, Q(H,H...) ~ (24/i 5 )M 2 /^ wh en /3 = -1/2. Hence, if M 2 > 0, 
the model can become singular when G — > _ (Type II singularity). 

• For N — 1 and n = 3, Q(H,H...) ~ -6102/63 when /3 = -3/7. Thus, if M2 < 0, the 
model can become singular when G — > 0~ (Type II singularity). 

In a certain sense, the model F\(G) in Eq. (I3.38P is a particular case of Eq. ( I3.39p . For 
large values of G, it tends to a constant with velocity being zero, so that it is impossible to 
find singularities (it is well known that R + constant is free of singularities, according to the 
ACDM model). Nevertheless, similarly to the above, a Type II singularity can occur when 
G -> 0- for n > 1. For example, if n = 2, Q(H, H) ~ (24/i 5 /6§)(&ia a - ai& 2 ) for /3 = -1/2. 
If 6ia 2 — a x 6 2 > 0, the model can become singular when G — > 0~. 
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With regard to F 3 (G) in Eq. (I3.40p . it is interesting to find the conditions on m, n, a 3 
and 63 for which we do not have any type of singularities. When G — > ±00 or G — > 0~, it 
is possible to write this model in the form F(G) ~ aG 1 , which we have investigated on in 
the preceding subsection. We do not consider the trivial case n — m. The no-singularity 
conditions follow directly from the results of the preceding subsection as complementary 
conditions to the singularity ones: 

• Case (A): n > 0, m > 0, n / 1 and m/1. We avoid any singularity if < n+m < 1/2 

and 0363 < 0; n + m > 1/2, n > 1 and a 3 > 0; n + m > 1/2, 2/3 < n < 1 and a 3 < 0; 
n + m > 1/2, < n < 2/3 and if n = 1/2, a 3 > 0. 

• Case (B): n > 0, m < and n / 1. We avoid any singularity if0<n<l/2 and 
03 < 0; n > 1/2, n + m > 1 and a 3 b 3 > 0; n > 1/2, 2/3 < n + m < 1 and 0363 < 0; 
n>l/2,n + m<2/3 and if n + m — 1/2, a 3 b 3 > 0. 

• Case (C): n < 0, m > and m 7^ 1. We avoid any singularity if m + n > 1/2; 
m + n < 1/2 and 0363 < 0. 

• Case (D): n < and m < 0. We avoid any singularity if a 3 < 0. 

We end this subsection considering the following realistic model, again for n > 0, 

F,(G)=G ^ n + b ' (3.46) 

Since for large G, one has F&(G) ~ ai/a 2 G a and for small G, one has F^G) ~ fei/^G", the 
analysis of Subsection III. B leads to the absence of any type of singularities for 

12 . 

- < a < - . (3.47) 

In fact, for this range of values, the asymptotic behavior of the right-hand side of Eq. (13. 8p 
is different from the asymptotic behavior of its left-hand side on the singularity solutions. 
Thus, Eq. (I3.46P provides an example of realistic model free of all possible singularities 
when Eq. (I3.47P is satisfied, independently of the coefficients. Moreover, this model suggests 
the universal scenario to cure finite-time future singularit y. Adding above model to any 
singular Dark Energy (in the same way as adding i? 2 -term 20, 
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results in combined 



non-singular model. Hence, unlike to convenient DE which may be singular or not, (non- 
singular) modified gravity may suggest the universal recipe to cure the finite-time future 
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singularity. In this respect, modified gravity seems to be more fundamental theory than 
convenient DEs. 



IV. FINITE-TIME FUTURE SINGULARITIES IN F{R, G)-GRAVITY 

In this section, we consider the finite-time future singularities in J-(R, G)-gravity. We 
reconstruct the ^(R, G)-gravity models with producing the finite-time future singularities. 
Furthermore, we examine the finite-time future singularities in a realistic model of ^(R, G)- 
gravity. 



A. Formalism 



We study the pure gravitational action of ^(R, G)-gravity, i.e., the action in Eq. (12.11) 
without /^matter- In this case, it follows from Eqs. ( 12. 6 P and (12. 7\i that the EOM of F(R, G)- 



gravity are given by 



35| 



24H 3 P G + m 2 p R + m P R + (J 7 - RP R - gp g ) = o , (4.1) 

8H 2 P G + 2P R + 4HP' R + 16HP' G (H + H 2 ) 

+P R {AH + QH 2 ) +F- RP R - GP G = . (4.2) 

In the case of pure gravity, these two equations are linearly dependent. 

Now, similarly to the previous section, by using proper functions P(t), Z(t) and Q(t) of a 
scalar field which is identified with the time t, we can rewrite the action in Eq. (12. ip without 

■^matter to 

S = i / ( P ^ R + Z ^ G + ■ ( 4 - 3 ) 

By the variation with respect to t, we obtain 

P'(t)R + Z'(t)G + Q'(t) = 0, (4.4) 

from which in principle it is possible to find t = t(R,G). Here, the prime denotes differen- 
tiation with respect to t. By substituting t = t(R,G) into Eq. ( 14. 3p . we find the action in 
terms of ^(R, G) 

J r (R,G)=P(t)R + Z(t)G + Q(t). (4.5) 
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By using the conservation law and Eq. (14. ip . we get the differential equation 

P"(t) + Ag 2 (t)Z"{t) - g(t)P'(t) + (8gg - Ag\t))Z'{t) + 2g(t)P(t) = , (4.6) 



where we have used the expression of the scale factor in Eq. (I3.14p and the Hubble parameter 
H(t) = g{t). By using Eq. (jPj) . Q(t) becomes 



Q(t) = -24g 3 (t)Z\t) - Qg 2 {t)P{t) - 6g(t)P'(t) . 



(4.7) 



In general, if P(t) ^ 0, J^(R, G) can be written in the following form: 



F(R,G)=Rg(R,G) + f(R,G) 



(4.8) 



where g(R, G) ^ and f(R, G) are generic functions of R and G. From Eqs. (14.11) and (14. 2p . 
we obtain 



Pcff 



2K 2 g(R, G) 



1 1 H*P G + 6H 2 [ R d -i^l + ^M)') + mP R 



dR 



dR 



+(T-RP R -GP C 



G) 



and 



PcS 



m 2 T G + 2P R + AHF' R + 16HP G (H + H 



2K 2 g(R,G) 

+ f R dg(R,G) , #(i?,G) 



+ 



(4ij + 6# 2 ) + J 7 - itF^j - GP G 



(4.9) 



(4.10) 



respectively, where p eS and p e g- are given by the expressions in (12.51) . As a consequence, we 
recover the same formalism of Sec. Ill as 



G(H,H...) = p cS - wp cS 
1 



2K?g(R,G) 



dg(R,G) | df(R,G) 



QH 2 (l + w) +4H 



dR dR 
+HP' R (A + 6w) + 8HP G \2H + H 2 {2 + 3w) 



+ 2p R + 8H 2 P G \, (4.11) 



where w is the constant EoS parameter of matter. The use of this equation requires that 
g(R, G) ytz on the solution. The equation for Q(H, H...) is given by Eq. (13. 8p . 
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B. Finite-time future singularities 

We examine the J-(R, G)-gravity models in which the finite-time future singularities could 
appear. 

1. Big Rip singularity 

First, we investigate the Big Rip singularity. If /3 = 1 in Eq. (12.81) with Hq = 0, we have 

H = T^-^ (4-12) 
6h 

R =u 7^( 2/i + 1 )' ( 4 - 13 ) 

(*o - ty 

24-h 3 

G ^M (1 + ,I) - (414) 

with h > 0. A simple (trivial) solution of Eq. (14.61) is given by 

P(t) = a{t - t) z , (4.15) 
Z(t) = 6(t -ty, (4.16) 

with a and 5 being constants, where x = 3 — h and z is given by 



l-h± Vh 2 -10h + l 
z± = W — . (4.17) 

Thus, the most general solution of P(t) is expressed as 

P{t) = «i(t - t) z+ + a 2 (t ~ t) z - , (4.18) 

where a% and a 2 are constants. From Eq. (14. 7p . we have 

^_ 2Ah 3 6(3-h) Qh ai (z + -h) Qha 2 (z^-h) 

m ~ (to - W 1 + (to-ty-^ + (to-*) 2 -*- ' ( } 

Under the condition < h < 5 — 2\/6 or h > 2 + a/6, the solution of J-(R, G) (by absorbing 
some factor into the constants) is given by 

F(R, G) = a x R x ~ z + 12 + a 2 R x ~ z ' t2 + 5G^ . (4.20) 

If 5 = 0, we find a well-known result of /(i?)-gravity. G^ is an invariant of the Big Rip 
solution in a F(G)-gravity theory and it is a solution in a general J-(R, G)-gravity theory. 
Note that 1 - z± ^ 1. 
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Another exact solution of Eq. (14. 6 j) is given by 

Pit) 



Z(t) 



(t - ty 

5 



(t - ty ' 

where 5 and x are constants, z = x + 2 and a is given by 

4h 2 5x(h -x-3) 
x 2 + (5 — + 6 



(4.21) 
(4.22) 



a 



(4.23) 



From Eq. (14.71) . we find 



Q(t) 



6h 



(to - ty +A 

The solution of Eq. ( 14. 4p is given by 



[Ah 2 8x + a(x + 2 + h)] . 



(4.24) 



t -t = f(R,G) 

-a(x + 2)R ± v/a 2 (x + 2) 2 R 2 + 24/i [4/i 2 5x + a(x + 2 + /i)] (x + 4)6xG]^ /2 



25xG 



(4.25) 



with i^O and 5^0. 

To have real solutions, we must require that the arguments of the roots in Eq. (I4.25P are 
positive. For h > 0, the principal cases are as follows: 



Case 
Case 
Case 
Case 
Case 
Case 
Case 
Case 



1) : x > 0, S > 0, 1 + x < h < x + 5 + f. We must use the sign + in (I4.25P 

2) : — | <£<0, 5 < 0, h > x + 1. We must use the sign +. 

3) : — 4 < x < — |, <5<0, h > x + 5 + -. We must use the sign +. 

4) : a;>0, <5<0, x + 5 + ^>h>l + x. We must use the sign — . 

5) : — | < x < 0, S > 0, h> x + 1. We must use the sign — . 

6) : —4 < x < — |, 5 > 0, h > x + 5 + -. We must use the sign — . 

7) : x = —4, 5 > 0. We must use the sign — . 

8) : x = —4, 5 < 0. We must use the sign +. 
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The solution of J-{R, G) is given by 

HR, G) = {nR ; Gyr , R + (7(^ G - {m %)^ ^ Sx + a(x + 2 + "W • < 426 > 

where f(R,G) is given by Eq. (I4.25p . This is an exact solution of EOM in Eqs. (14. ip and 
(14. 2 p for the Big Rip case. 

We show several examples. In the case a = 1 and x = —2, we find 

m ,G) = R+^y^SvG, Ml, (4.27) 
(1 - h) 

which is in agreement with the result of the previous section. 

If a = and x = h — 3 (this case corresponds to the cases (l)-(6) presented above), we 
find 

F{R, G) = SG^ , 5^0, (4.28) 

which is equivalent to Eq. (14.201) with ot\ = a 2 = 0. 
If x = —4, the result is given by 

16h 4 8 

G 



F(R,G) 



5^0. (4.29) 



l + 2/i 2 ) 2 

Hence, if for large values of R and G, ^(R, G) ~ ±a =F 5(R 2 /G) with a > and 5 > 0, the 
Big Rip singularity could appear. 

If x = h — 1, the solution becomes (by absorbing some constant) 

i-fe 
R\ ~ 



F(R,G) = 6G^-J , 5^0, Ml. (4.30) 

Thus, if for large values of R and G, F(R, G) ~ 5G"< / R^ 1 with 5 ^ and 1/2 < 7 < 1 or 
1 < 7 < +00, the Big Rip singularity could appear. 
Furthermore, it is possible to verify that the model: 

F(R,G)= 1 — 1 (4.31) 

with 7 being a generic constant, is a solution of Eqs. ( 14. ip and (14. 2 p in the case of the Big 
Rip singularity (/3 = 1) for some value of h. In general, it is possible to obtain solutions 
for h > if m > 0, n > and m > n. For example, the case n = 2 and m = 3 realizes 
the singularity in h = 5; the case n = 1 and m = 3 realizes the singularity in /i = 4 + \/19 
and so forth. This is a generalization of Eq. ( I4.30p . Note that we do not recover a physical 
solution for m — — 1 and n = —2 because in this case h = —3. For a similar kind of model 
F(R 2 /G) which produces the Big Rip singularity, see Eq. (I4.29p . For m = or n = 0, we 
recover Eq. (I4.20p . 
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2. Other types of singularities 



Next, we study the other types of singularities. We consider the case in which H is given 



by 



H 



(t - ty ' 



An exact solution of Eq. (14. 6 p is given by 



P(t) = -\(4h 2 )(t -t), 

z(t) = \(t -t) 2 e+\ 



(4.32) 

(4.33) 
(4.34) 



where A is a generic constant. The form of Q(t) is given by 

Q(t) 



24h A \ 48h 3 f3 
+ 



(4.35) 



(f -f)a/>-i (t -ty 

For = 1, we find a special case of Eq. ( I4.26p . For (3 > 1, we obtain the asymptotic real 
solution of Eq. P~4l) : 

1/2/3 

(4.36) 



t -t = f(R,G) = 2 1 ^ 
The form of F(R, G) is expressed as 



h 2 R + ^/h 4 R 2 + 6/i 4 (4/3 2 - 1)G 



[l + 2/3)G 



F{R, G) = -Ah 2 X(f(R, G))R + \(f(R, G) 1+W )G + 24h 4 \(f(R, G) 



\1-2B\ 



0>1. (4.37) 



This is an asymptotic solution of Eq. (13. 8p when (for (3 > 1) 



1 



K 



2H + 3(1 + w)H 2 



3(1 +w)h 2 



K 



(t - t) 



-23 



(4.38) 



In the case /3 > 1, the form of J-(R, G) is written as 

aG 



F(R,G) ~ A 



R 



a > 0, 7>0, A ^ 0. 



(4.39) 



This is the asymptotic behavior of a J-(R, G) model in which a "strong" Type I singularity 
(P > 1) could appear. By taking g(R, G) = ~fG m /R n+1 and using Eqs. gS]) and (l4TTU]l . it 
is possible to verify that for the model 

G' 



F(R,G)= 1 



R n 



(4.40) 
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the function Q(H,H...) in Eq. (14. lip is given by 

g(HH ^ JhH2m- „-!)(! + W ) 

which is, under the condition 2m — n — 1 > 0, an asymptotic solution of Eq. (13. 8 P in the case 
of f3 > 1. Thus, in the model J-"(-R, G) ~ jG m /R n with m > (n + l)/2 the Type I singularity 
could appear. This point has important consequences because it is possible to see that the 
theories F(R) = R n with n > 1 or F(G) = G m with m > 1/2 can become singular. 

To find other models, we can consider the results of Sec. III. The Type I singularities 
correspond to the asymptotic limits for R and G 

R ~ 12H 2 , G ~ 2AH A . (4.42) 

These are two functions of the Hubble parameter only, so that 

lim 24 (—} = lim G . (4.43) 

M-to \ 12 / t-+t 

If we substitute G for R in Eq. (13.301) by taking into account Eq. ( 14.431) . we obtain a zero 
function (this is because Eq. ( 13.30}) is zero on the singularity solution). If we substitute G 
for G/R, however, we obtain the following model: 

F{R, G) = R - ^ . (4.44) 

This is an asymptotic solution of Eq. ( 13. 8 P such as Eq. (I4.38p . Thus, there appears Type I 
singularity for J-"(-R, G) ~ R — a(G / R) with a > 0. 

In the case of H = h/(to — t) 13 with < 1, it is not possible to write G and R like 
functions of the same variable (H or the same combination of H and H). Nevertheless, if 
we examine the asymptotic behavior of G and R, we have 

6hB 

R ~ (4 45) 

24h 3 f3 

and 



G - u *s P+ i . ( 4 - 46 ) 



G 2/3 

_ ~ G 3^+i (4.47) 

If we use G/i? for G in Eq. ( 13.341) as in Eq. (14.4Tj) . we see that the asymptotic time dependence 
in Eq. (13.81) for (3 < 1 is the same: 



1 

K 2 



2H + 3(1 + w)H' 



a 7 
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+ ■ (4-48) 



Under this consideration, it is possible to derive a J-(R, G)-gravity theory (by setting some 
parameters) from Eq. f)3.34p as 



JP(R,G) = R+^, (4.49) 

in which the other types of singularities appear. Thus, in this model (^(R, G) ~ R+a(G/R) 
with a > 0) the Type II, III and IV singularities could appear. Then, by substituting G for 
R we get 

(1 + 0) 2/3 

?{R,G) -R-5 y _ i- \R\m , 5>0. (4.50) 

KP / 

This is a well-know result. In the model J-"(i? — > oo) ~ R + aR 1 , for < 7 < 1 and a > 0, a 
Type III singularity could appear. In the model JF(R — > —00) ~ -R + a|-R| 7 , for —00 < 7 < 
and a > 0, a Type II singularity could appear. In the model F(R — > 0~) ~ R + a|i?| 7 , 
for 2 < 7 < +00 (7 7^ 2n/(n — 1), where is a natural number) and a < 0, a Type IV 
singularity could appear. (In the Big Rip case, we have found exact solutions. This kind of 
reasoning is therefore inapplicable.) 

C. Realistic model of T(R, G)-gravity 

We study the following realistic model of T{R, G)-gravity: 

J"(-R, G) = a\G n + a 2 R m + - ° 5 - , (4.51) 

where diii = 1, • • • ,5) are constants and n(> 0) and m(> 0) are positive constants. For 
large values of R and G, we have 

7{R,G) ^ ai G n + a 2 R m . (4.52) 



In the specific case n > 3 and m = (1/2) (1 + 2n + a/3 — 12n + 4n 2 ) (in which, m > 
(7 + V3) /2), from Eq. (jODjl we see that the Big Rip singularity could occur. To find other 
singularity solutions, we investigate the asymptotic form of Q(H,H...) in Eq. (14. lip and 
require the consistence with Eq. (13.81) of Eq. ( 14. 81) . The behavior of Eq. (14.111) takes two 
different asymptotic forms which depend on the parameter of as follows: 



22 



Case of j3 > 1: In the limit t — > to, we find 



Q(H, H...) 



1 



g(R,G) 



a 



V 

■J R 



+ 



+7 



+5- 



dg(R,G) df(R,G) 



(to - ty? 



+ e- 



dR 
V 

•J r 



G 



dR J (t - i) 2/3 



+ (P R + V 



G 



(4.53) 



[to-ty (t -tye '(t -tye 

where e and rj are constants. To realize a I Type singularity, from Eq. ( 13. 8 j) we must 
have 

3(l + w)/i 2 



G(H,H.. 



(4.54) 



/? 2 (t - 1) 2/3 ' 

Hence, if for G ~ l/(t - t) 4/3 and R ~ l/(t - t) 2/3 with (3 > 1, the highest term of 
Eq. (I4.53P is proportional to l/(to — £) 2/3 , it is possible to have a Type I singularity. 
As in F(G)-gravity, this condition is necessary and not sufficient. 



Case of /3 < 1: In the limit t — > to, we obtain 
1 



Q(H,H. 



g(R,G) 



a 



X R 



(to ~ t) 



+ 



X G 



(to - 1) 



3/3+1 



+ 7 



dg(R,G) 



(to ~ t) 



df(R,G) 
dR 



V 

+ 6- 



+ e- 



^G 



CPr + v-, T ' G 



(to-ty +1 (to-ty (to-ty^ * '(t -tye 

To realize this kind of singularities, from Eq. (13. 8p we must have 



dR 

. (4.55) 
(4.56) 



K\t -ty^- 

Thus, if for G ~ l/(t - t) w+l and R ~ l/(t - t) /3+1 with /3 < 1, the highest term of 
Eq. (I4.55P is proportional to 1/ (t — ty +1 , it is possible to have a Type II, III or IV 
singularity. Also this condition is necessary and not sufficient. 



For the model in Eq. (I4.5ip . if m and n are positive numbers (F(R — > oo,G — > oo) ~ 
a\G n + a 2 -R m ) and m = 2n, the asymptotic behavior of Eq. (I4.53P (in this case g(R, G) = 
a 2 R m - r ) when R ~ 12h 2 /(t - t) 2/3 and G ~ 24/i 4 /(t - t) 4/3 is proportional to l/(i - t) 2f3 
and therefore it is possible to realize the Type I singularity. As a consequence, we get 



Q(H, H...) 



3(1 + w)h 2 



m 



1 + 



(m — 2)eii 
a 2 



(4.57) 



/<t -*) 2/3 

To have the consistence with Eq. (I4.54|) . we find that if 1 < m < 2 and 01/02 < 0, the Type 
I singularity could appear (for example, if m = 1, n = 1/2, 02 = 1 and ai < 0, we recover 
the case of Eq. (|3.30|) ). The same result is obtained if m > 2 and a 1 /a 2 > 0. 
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We see that for m > 0, n > and — 1 < /3 < 1, j 7 ' R /g(R,G) ~ (t — t) 2 (like above, 
g(R,G) = a 2 R m ~ 1 ) and Q{H,H...) in Eq. (jOgD diverges faster than (t -t) - " -1 , so that in 
order to find Type II or III singularities, we must have J 7 ^ = 0. In general, this is true if 
m = 1 and we can recover the results in Sec. Ill for F(G) -gravity. 

When R — > 0~ and G — > 0~, for m > and n > 0, J-(R, G) behaves as 

T{R -» 0", G -» 0") ~ ^ ° 5 „ ■ (4.58) 

In this case, if /3 < -1, Eq. (#(i?, G) = a 5 /(a 3 G n J R + a 4j R m+1 )) diverges and Eq. (jOSjl 

becomes inconsistent, so that the model is free of Type IV singularities. 



V. CURING THE FINITE-TIME FUTURE SINGULARITIES 

In this section, we discuss a possible way to cure the finite-time future singularities in 
_F(G)-gravity and J-(R, G) -gravity. In the limit of large curvature, the quantum effects 
become important and lead to higher-order curvature corrections. It is therefore interesting 
to resolve the finite-time future singularities with some power function of G or R. 

A. i ? (G)-gravity 

First, we consider F(G)-gravity. If any singularity occurs, Eq. ( 13.81) behaves as 

'_ w^)# + m {to _ t) -2 BigRip 

G(H,H...)~ I _3(i±^(t _t)-2/3 /3>1 (Type I) (5.1) 
k -^r(*o - t)-f*-i $ < 1 (Types II, III, IV ) 

The singularities appear in two cases: G — > ±oo or G — > 0~. 

(i) Case of G — > ±oo 

Suppose that for large values of G, 

R + F(G^±oo) — >R + ^G m , m^l, (5.2) 

with 7 7^ 0. One way to prevent a singularity appearing could be that the function Q(H, H...) 
becomes inconsistent with the behavior of Eq. (15.11) . In general, Q(H,H...) must tend to 
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infinity faster than Eq. (15.11) . For H = h / (to — t) (Big Rip), we have 

Hence, if m > 1/2, we avoid the singularity. Nevertheless, there is one specific case in which 
the Big Rip singularity could occur. If m — (1 + h)/A, Q(H,H...) is exactly equal to zero, 
so that (for example) the following specific model admits the Big Rip singularity: 

This is because the power function G m is an invariant with respect to the Big Rip singularity 
of G 1 / 2 . If for large values of G, F(G) ~ «G 1/2 , we can eliminate the Big Rip singularity 
with a power function / yG m (m > 2) only if a > 0. 

For H = h/(to — t)P with (3 > 1 (Type I) and the behavior in Eq. ( 15. 2ft . we find 



Also in this case, if m > 1/2, we avoid the singularity. For example, R + F(G) = R + ayG + 
'jG 2 with a > is free of Type I singularities, while if a < 0, the Big Rip singularity could 
appear. 

For H — h/(to — t)P with < (3 < 1 (Type III) and the behavior in Eq. (15. 2p . we obtain 

6( g -g-)^ (to _ t) ^ +1)+(1 - a - < 5 - 6 > 

If m > 2/3/(3/3 + 1) (i.e. m > 1/2), we avoid the singularity. 

Also for H = h/(to — t) 13 with — 1/3 < /3 < (Type II, G — > — oo), we have to require 
the same condition. For example, R + a\G\^ + 7G 2 with £ < 1/2 is free of Type I, II and 
III singularities. 

(ii) Case of G ->■ ID- 
Suppose that for small values of G, 

R + F(G -> 0~) — >■ R + 7G m , (5.7) 

with 7^0 and m being an integer. For H — h/(to — t) 13 with /3 < —1/3 (Type II and IV 
singularities), we get 
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which diverges and hence becomes inconsistent with Eq. (15. II) if m < 2/3. For example, 
a:|Crp + lG~ l with Q > 2/3 is free of Types IV singularities. 

As a result, the term ^G m with m > 1/2 and m ^ 1 cure the singularities occurring 
when G — > ±oo. Moreover, the term '-fG" 1 with m < and m being an integer cure the 
singularities occurring when G — > 0~. 

In /(i?)-gravity, by using the term ^R" 1 , the same consequences are found. The term 
jR" 1 with m > 1 cures the Type II and III singularities. On the other hand, the term jR 171 
with m < 2 cures the Type IV singularity. 

Note that ^G" 1 or jR" 1 are invariants with respect to the Big Rip solution (see Eq. f)4.20p ). 
so it is necessary to pay attention to the whole form of the F(G) or /(i^)-gravity (see 

Eq. 453D). 

It is also possible to cure the singularities in a F(G)-gravity theory with the power 
functions of R and a /(-R)-gravity theory with the power functions of G. To do it, it is 
useful to take into account that G diverges as R 2 in the Type I singularity solutions, at least 
as R in the Type II, and as R 2 in the Type III, and G tends to zero at least as R 3 in the 
Type IV (proved by the fact that when f3 < 1, G ~ R 3/3+1 ^ +1 ). We show several examples. 

• R + G n + R m ~ R + R m on the asymptotic limit of the Types I and III singularity 
solutions (if they exist) when m > 2n. 

• R + G n + R m ~ R + R m on the Types II singularity solutions if m > n (on this kind 
of solutions R tends always to infinity, while in some cases, for —1 < < —1/3, G 
tends to zero). 

• l/(G n + R m ) ~ 1/R m on the Types IV singularity solutions (for which R — > 0~ and 
G (T) if m < 3n. 

Thus, the singularity solutions found in Sec. Ill for F(G)-gravity can be cured by the term 
^R m with m > 1 for Type I, II and III singularities and that with m < 1 for Type IV 
singularity. 

We mention the Type I singularities with /3 > 1. We have shown that the model R + ~fG m 
with m > 1/2 (m ^ 1) is free of Type I singularity because Eq. (13. 8p becomes inconsistent. 
Nevertheless, it follows from Eq. f)4.40p that the models G m with m > 1/2 and R n with 
n > 1 can show the Type I singularity in the asymptotic limit. This means that when t is 
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very close to to, the term G m or R n is dominant over R, R + G m ~ G m (or R + R n ~ R n ) 
and therefore the Type I singularity could appear. Hence, the important point is whether 
the model can approach to very large values of R or G with a non-singular metric (which is 
not admitted) and then become singular because R is negligible. It depends on the form and 
the dynamics of the model and the value of m or n. If m ^> 1/2 or n ^> 1, the singularity 
could appear more easily. Thus, in order to avoid the singularity solutions, it is better to 
choose m and n as m > 1/2 and as n > 1, respectively, but m and n are not very large. 
This does not hold in the other types of singularities. The theory R + G m with m < and 
m being an integer (or R + R n with n < 2) is free of Type II, III and IV singularities as the 
theories G m or R n . 



B. F{R, G)-gravity 

Next, we study ^(R, G)-gravity. In the general F(R, G)-gravity, in order to avoid the 
singularities with power functions, we must require that the EOM (14. ip and (14. 2[) are incon- 
sistent on the singularities solutions. Within the framework of J-"(R, G)-gravity, we can use 
the terms such as G m /R n to cure the singularities. The singularities appear in the following 
three cases: (a) R — > ±oo, G — > ±oo (Types I, II, III), (b) R — > — oo, G — > CT (Type II for 
-1 < < -1/3), and (c) R 0~, G (T (Type IV). 

We investigate general cases. Suppose that for large values of G and R, 

G m 

J^(R^oo,G^oo)^R + 1 —, (5.9) 

with 7^0. In the case of the Big Rip singularity, in which H is given by Eq. (I3.17p . 
Q(H,H...) in Eq. (14. lip diverges as 



G(H,H...)~ {tQ _^ m . 2n - (5-10) 

Thus, if m > (n + l)/2, we avoid the singularity. Nevertheless, there is the possibility that 
Q(H,H...) is exactly equal to zero and the Big Rip singularity could occur (see Eqs. (I4.30p 
and (I4.3ip in the case of m — n + 1). Hence, the whole form of ^(R, G) as well as its form 
in the asymptotic limit must be examined. 

In the case of Eq. (I3.27P (Type I), Q(H,H...) diverges as 

OL 

G(H, H...) ~ — — . (5.11) 
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Also in this case, if m > (n + l)/2, we avoid the singularity. Similarly to the above, however, 
if m 3> 1 and n C 1, the asymptotic limit of T{R, G) in Eq. (15. 9p behaves as / yG m /R n and 
therefore the Type I singularity could occur (see Eq. (I4.40p ). 

As a consequence, we can avoid the Type I singularities if the asymptotic behavior of the 
model is given by Eq. (15.91) and its asymptotic form has the power functions 

G m R n , (5.12) 



or 

G m n + 1 

~R™ 



m>— — , (5.13) 



with m and n being positive integers. 

Now, suppose that when H = h/(t Q — tf with /3 < 1, the asymptotic limit of F{R, G) 
becomes 

r(R,G)-+-y— . (5.14) 

For < 1, Q(H,H...) behaves as 

which diverges faster than (to — t) ~P~ l and therefore the Type II, III and IV singularities are 
always avoided for any value of m and n. The same scenario to cure the future singularity 
by adding the non-singular modified gravity maybe applied here again. 

VI. CONCLUSION 

In the present paper, we have investigated the finite-time future singularities in F(G)- 
gravity and T{R, G)-gravity. We have reconstructed the F(G)-gravity and ^(R, G)-gravity 
models in which the finite-time future singularities may occur. It has been demonstrated 
that all four types of finite-time future singularity may emerge for a variety of the above 
models with the effective quintessence/phantom EoS behavior in the same qualitative way 



as for convenient DEs where also all four types of future singularity may occur [17H19]. This 
provides the explicit demonstration that whatever the effective DE model (convenient one or 
modified gravity) is, it may lead to singular future universe. Moreover, the future singularity 
may manifest itself as radius singularity for spherically-symmetric spaces. This may cause 
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instabilities for black holes [18| and relativistic stars |27l. l29l. l3ll. 39 |. Other imprints of the 
singular future universe to current cosmology may be searched as well. 

However, there exists fundamental qualitative difference between convenient DE and 
modified gravity. It turns out that sometimes it is possible to solve the singularity issue 
taking account of quantum gravity effects (see Big Rip singularity resolution in Ref. 40]) 
or by the coupling of DE with Dark Matter (DM) (some fine-tuning of initial conditions 



may help to resolve Type II or Type IV future singularity [41]). Nevertheless, quantum 
gravity account is effectively the modification of gravity. Moreover, it is only modified 
gravity (actually, its additional modification as we have demonstrated on the example of 
non-singular F(G)-model in Subsection III. C) may suggest the universal scenario to cure 
any finite-time future singularity. This is achieved by adding such non-singular theory to any 
DE containing future singularity in its evolution. Furthermore, such additional modification 
may always be made by terms which are relevant only in the early universe and are typical 
as quantum gravity corrections. Hence, modified gravity suggests the universal scenario to 
protect the future universe from singularity while not destroying the attractive cosmological 
properties of specific DE alternative gravity like its viability with local/cosmological tests if 
exists. This may be considered as powerful theoretical argument in favor of the consideration 
of such theories as DEs. 
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Appendix A: Simple model of /(i?)-gravity 

There is a great diffusion of modified gravity models which for large values of curvature 
tend to a constant and imitate the ACDM model. Nevertheless, the existence of derivatives 
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of the modified function, which is very small but different from zero, involves the possibility 
of singularities for < 1. In this Appendix, we explore the finite-time future singularities 
in a simple model of /(-R)-gravity. The action of /(i?) -gravity is given by 

1 

tq. (ED) with F(R,G) = R + f(R). We examine the 



S — d A xy/^g 



, (R + f(R))+£ 



(Al) 



which corresponds to the action in 
following simple /(i?)-gravity model 



42] which reproduces the current accelerated expansion 



of the universe and imitate a cosmological constant for large values of curvature: 

f(R) = a(e~ bR - 1) , (A2) 

where b (> 0) is a positive constant. In what follows, we consider the pure gravitational 
action of /(i2)-gravity in Eq. (1A1I) without £ m atter- 



The finite-time future singularities in /(-R)-gravity have been discussed in Ref. 2l| . 
from which we propose the asymptotic expression of Q(H,H,...) when < 1, similarly 
to Eq. (ET351) : 



Here, the prime denotes differentiation with respect to R. When /3 < 0, the asymptotic 
behavior of R is given by 

R ~ Qhf3{t - t)-?- 1 . (A4) 
We assume R > 0, so in this case h must be negative. If H behaves as 

H ~ h(t — t)~~P + H , (A5) 

as Eq. ( 12. 8p . H can be still positive in the limit t — > to. For large values of R (— 1 < /3 < 1), 
f(R) in Eq. tend to -a and 

f'(R)~-abe- bR , f"(R) ~ab 2 e- bR , f"'(R)~-ab 3 e~ bR . (A6) 

To obtain a singularity, the highest term in Eq. (1A3I) must be divergent as l/(£o — In 
order to check it, it is convenient to develop the exponential function in power-series. The 
third term of Eq. (1A3|) behaves as 

" f"(R) - ~ u ■ (A7) 
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For n = 2/(/3 + 1), the highest term of the denominator behaves as (to — If j5 — > — 1 + , 
n — > oo and this is just the asymptotic value of n of the highest term in Eq. (1A7I) . A 
similar argument is valid for the last term of Eq. (IA3[) . whereas the first and second terms of 
Eq. (IA3p tend to a constant and zero, respectively. Thus, a Type II singularity could occur in 
the model in Eq. (1A2I) for large values of R. Note that /(i?)-gravity unifying the early-time 
inflation with late-time acceleration as proposed in Ref. |43| turns out to be non-singular 
due to the presence of R 2 term. 



Appendix B: Asymptotic behavior of singular models 



In this Appendix, we discuss the asymptotic behavior of singular models. In Sec. IV A, 
we have shown that in principle it is possible to write a general F(R, G)-gravity theory in 
the form in Eq. (14.81) . In this case, from Eqs. (13. 9p and (14.1 ip we obtain 

1 



G(h,h...) = - 



2H + 3(1 + w)H' 



(Bl) 



where 



Q{H,H...) 



1 



2K 2 g(R,G) 



'l+w)(F-RT R -GF' c 



G) 



{R dg(R,G) | df(R,G) 



d R ' d R J^(1 + W ) + 4H 
+HP R (A + 6w) + %HP G \2H + H 2 {2 + 3w) 



2J-' R + SH 2 J-q 



(B2) 



It is clear that in the case of F(R, G) = R+F(G), by taking g(R, G) = 1 and f(R, G) = F(G) 
in Eq. flUD, Eqs. Q£5D and (l3~T0|) in Sec. Ill A is recovered. 

To verify the existence of singularities, it is very useful to control the consistence of 
Eq. (IBip . On the exact solutions (this is the case of the Big Rip solutions), this check 
is independent of the choice of g(R,G) in Eq. (|4.8|) . We must carefully verify only that 
g(R, G) 7^ on the singularity solution (or -equivalently- if we use this equation to find 
singularity solutions, it is possible to lost some solutions in which g(R,G) = 0). When we 
check asymptotic solutions, a problem could appear. It is based on the confrontation between 
the asymptotic behaviors of the right-hand and left-hand sides of Eq. ( IBip . Furthermore, in 
this case there is a problem if all the terms of F(R, G) have the same asymptotic behavior 
on the singularity solution. For example, in order to verify the existence of the Type I 
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singularity solution in the model R — aG/R with a > (see Eq. f)4.44p ). it is indifferent to 
take g(R,G) = 1 and f(R,G) = -aG/R or g(R, G) = -aG/R 2 and f(R, G) = R because 
on the singularity solution the asymptotic behavior of R (R ~ (to — t) 2 ) is the same as 
G/R. Nevertheless, if the terms do not have the same asymptotic behavior (for example, 
J-(R, G) = R + R 2 + R 3 ...), it is necessary to be careful in the choice of g(R, G) to substitute 
into Eq. (1B2I) . The mechanism is the following: The right-hand side of Eq. ( 1B1|) behaves 
as R, while the left-hand side of Eq. (1B1[) is proportional to l/g(R,G). Automatically, in 
the asymptotic limit all the terms smaller than Rg(R, G) (and also their derivatives) are 
neglected. As a consequence, when we use the EOM in the asymptotic limit, we must choose 
g(R,G) as the coefficient of the smallest term which we want to consider. This is easy to 
do when the terms are completely different in the limit. For example, R 2 + 1 / R ~ R 2 when 
R — > oo (this is the principle that we have used in the study of the realistic models shown 
in the present paper). The question is trickier when the terms of J-"(R, G) tend together to 
infinity or to zero with different velocities. In this case, the choice of g(R,G) depends on 
our target, if we want to verify the EOM in more or less strong limit. For example, let us 
consider the model ^(R, G) = R + R 2 . If we choose g(R, G) = 1 and f(R, G) = R 2 , we find 
that Eq. (IBlj) is inconsistent on the Type I singularity solution, so we can say that the model 
is free of this kind of singularity. Nevertheless, the choice of g(R, G) = R and f(R, G) = R 
is equivalent to neglecting the first term of F(R, G), so we are considering F(R, G) ~ R 2 
(strong limit when R — > oo). In this case, we find that the model could be affected by the 
Type I singularity. The physical meaning has been discussed in Sec. V. 
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